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ment of strangeness observed in S + Au collisions, interpreted previously as 
possible evidence for quark-gluon plasma formation, has its origins in non- 
equilibrium dynamics of few nucleon systems. A factor of two enhancement 
of A at mid-rapidity is indicated by recent pS data, where on the average 
one projectile nucleon interacts with only two target nucleons. There ap- 
pears to be another factor of two enhancement in the light ion reaction SS 
relative to pS, when on the average only two projectile nucleons interact 
with two target ones. 

25.75+r ; 24.85. +r ; 24.85. +p ; 25.40.-h ; 25.40.Ve ; 24.10Jv 



I. INTRODUCTION 

The search for new states of dense nuclear matter is one of the most active areas of 
research in nuclear physics JlJ , [0 . Enhanced strangeness production in ultra- relativistic 
heavy ion collisions was suggested long ago 0] as a signal for quark-gluon plasma formation 
H, ||, and has been observed at both the AGS and SPS. There is extensive data from 
both the SPS at CERN fl7J- and the AGS at BNL PJ- p§ on strangeness yields from 



reactions ranging from elementary p + p to p + At and As + At for targets ranging up 
to At ~ 200 and beams up to = 30. Detailed rapidity and transverse momentum 
spectra of (K + ,K~ , K®,A ,A) are available and spectra of S~ and even Q~ are becoming 
available. In all cases their yield relative to pions or negative hadrons are larger in 
nucleus-nucleus than expected from geometrically scaled proton-proton collisions. New 
experiments with truly heavy ion projectiles are in progress with Au beams at BNL 
47j and with Pb beams at CERN (Pb(170AGeV) + Pb) |] and will soon extend 



considerably the data base. 

These and other data on nuclear reactions have stimulated the development of many 
hadronic transport models to address the problem of multi-particle production in nuclear 
collisions. These include Dual Partons Models(DPM) |48|| - |]54f , Quark Gluon String 
Models(QGSM) |5§- VENUS FRITIOF [§TJ- ATTILA [§§, HIJING @]- 



RQMD models |9|- 0, Parton String Model (PSM) [[T|,HIJET ]7|, Parton Cascade 



Model(PCM) |77| - |J75| . An excellent review and detailed comparison of the models is given 
by Werner in Ref. . 

At present no conventional explanation of the large enhancement of hyperons or anti- 
hyperons has been found. The Pomeron exchange picture has motivated the development 
of many of the above models with the Pomeron modeled in terms of colored strings. 
However, the string picture itself suggests the possibility of new dynamical mechanisms 
ranging from string fusion to color rope formation. Some of the above transport models 
like RQMD |74[ and VENUS |)0[ include such non-conventional mechanisms as default 



options. These proposed novel non- equilibrium dynamical mechanisms were shown to be 
able to reproduce many features of the observed strangeness enhancement |80| - f82fl , |72|1 , 



73] , [ |60|j . On the other hand, there have been many attempts (see, e.g., review by Heinz 
in m p. 205c and references therein) to attribute the strangeness enhancement to the 
formation of an equilibrated fireball containing a quark-gluon plasma state (51,0. 

Therefore it appears that either non-conventional multi-particle mechanisms or the 
existence of a new form of matter seems to be indicated by the observed strangeness 
enhancement. Either case is of basic interest. The goal of the present study is to clarify 
which of these alternatives is more compelling. We use the HIJING model to 
perform a linear extrapolation of strangeness production dynamics from pp to AA taking 
into account essential nuclear geometry and kinematical constraints. At higher collider 
energies it includes pQCD semi-hard processes, but in the SPS range it reduces essentially 
to a hybrid version of the FRITIOF and DPM models. We use the VENUS model |Kj 



to estimate possible effects of final state cascading and new mechanisms of strangeness 
production in few nucleon processes. The non-conventional mechanism in VENUS4.13 is 
the occurrence of "double strings" which may form when one projectile nucleon interacts 
with two or more target nucleons. A double string is defined as a color singlet baryon 
configuration consisting of one projectile quark connected to two different valence quarks 
in the target via a three gluon vertex. In earlier versions of the model the parameterization 
of the vertex kinematics led to anomalously large baryon stopping power. In the present 



version, the double string phenomenology is constrained to reproduce the pA — > pX data. 



However, the new feature, see eq. 15.52 in ref. |S0" |, is the assumption that the probability 
for hyperon production in the fragmentation regions is enhanced by a factor of two relative 
to the single string rates. This enhanced strangeness production mechanism due to double 



strings is similar to that postulated in the color rope model |[3| and incorporated into 
the RQMD model. The hyperon enhancement in VENUS is however more confined to 
the fragmentation regions. 

Both HIJING and VENUS have been compared to a wide variety of data in pp , pA and 
AA collisions [66[] , [ 67j , |6(| . However, no systematic study of strangeness production at 
SPS-CERN energies were performed up to now. In addition, there have been substantial 
changes in the final published data relative to earlier comparisons to preliminary 
data 0, ||. In this paper, we calculate the rapidity and transverse momentum spectra of 
strange particles for pp , minimum bias collisions of pS,pAg,pAu and central collisions of 
S + S, Ag, Au, W at the energy of 200 AGeV and Pb + Pb at the energy of 170 AGeV. 
We focus special emphasis on the comparison with the data on pp, pS, and SS from 
Alber et al. |18| . That comparison reveals that much of the enhancement of strangeness 
in heavy ion collisions can be traced back to the enhancement of strangeness in the 
lightest nontrivial ion collisions, p + S. Our main conclusion based on these data is that 
the enhancement of strangeness observed in S + Au is therefore most likely due to new 
non-equilibrium multi-particle production mechanisms in processes involving few nucleon 
systems. 

This paper is organized as follows: A brief description of the HIJING Monte Carlo 
model and theoretical background are given in Section II. For a detailed discussion of 
the VENUS model, we refer to the review in ||60|| . In Section III, detailed numerical 
results with HIJING and VENUS for pp,pA and AA reactions at CERN-SPS energies 
( yfs ~ 20A GeV ) for strangeness production are compared to experimental data and 
other model predictions . Section IV concludes with a summary and discussion of results. 
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II. OUTLINE OF THE HIJING MODEL 



A detailed discussion of the HIJING Monte Carlo model was reported in references 
H- [H]. The formulation of HIJING was guided by the LUND-FRITIOF and Dual 



Parton Model(DPM) phenomenology for soft nucleus-nucleus reactions at intermediate 
energies < 20 GeV) and implementation of perturbative QCD(PQCD) processes in 
the PHYTHIA model |84j] for hadronic interactions. We give in this section a brief review 



of the aspect of the model relevant to hadronic interaction: 

1 . Exact diffuse nuclear geometry is used to calculate the impact parameter dependence 



of the number of inelastic processes [63 



2. Soft beam jets are modeled by quark-diquark strings with gluon kinks along the 
lines of the DPM and FRITIOF models. Multiple low p? exchanges among the end 
point constituents are included. 

3. The model includes multiple mini-jet production with initial and final state radiation 
along the lines of the PYTHIA model and with cross sections calculated within the 
eikonal formalism. 



4. Hadronization is performed via the JETSET7.2 algorithm |8J] that summarizes data 
on e + e~. 

5. HIJING does not incorporate any mechanism for final state interactions among low 
Pt produced particles nor does it have color rope formation. 

The rate of multiple mini-jet production in HIJING is constrained by the cross sections 



in nucleon-nucleon collision. Within an eikonal formalism |85j the total elastic cross 
sections a e i, total inelastic cross sections a in and total cross sections a tot can be expressed 
as: 

a el = n dfe 2 (l-exp(- X (M)) 2 (1) 
Jo 

roo 

a in = 7c db\l-exp{-2 X {b,s))) (2) 



o 



oo 
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a tot = 2 7T / db\l - exp(- X (&, s))) (3) 



o 



Strong interactions involved in hadronic collisions can be generally divided into two cat- 
egories depending on the scale of momentum transfer q 2 of the processes. If q 2 < Aq CD 
the collisions are nonperturbative and considered soft and modeled by beam jet frag- 
mentation via the string model. If q 2 ^> Aq CD the subprocesses on the parton level are 
considered hard and calculated via pQCD [B7 . 



In the limit that the real part of the scattering amplitude is small and the eikonal 
function s) is real, the factor 

0(M) = l-exp(-2x(M)) (4) 

can be interpreted in terms of semi-classical probabilistic model as the probability for an 
inelastic event of nucleon-nucleon collisions at impact parameter b which may be caused 
by hard, semi-hard or soft parton interactions. 

To calculate the probability of multiple mini-jet, the main dynamical assumption is 
that they are independent. This holds as long as their average number is not too large as 
is the case below LHC energies []S7f| . When shadowing can be neglected, the probability of 
no jets and j independent jet production in an inelastic event at impact parameter b, can 
be written as : 

g (b, s) = (1 - exp(-2 x,(b, s))) exp(-2 Xh(b, s)) (5) 

gj (b,s)= [2Xh ^ s)]J .exp(-2 Xh (b,s)) j>l (6) 

where Xs{b, s) -is the eikonal function for soft interaction, 2 Xh{b, s) -is the average num- 
ber of hard parton interactions at a given impact parameter, exp(— 2 Xs{b, s)) -is the 
probability for no soft interaction. Summing eqs. (5) and (6) over all values of j leads to : 

oo 

J2 g j (b,s) = l-e X p{-2x.(b,s)-2 X h(b,s)) (7) 

3=0 

Comparing with eq.(HI) one has : 
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X(b, s) = Xs{b, s) + Xh{b, s) (8) 

Assuming that the parton distribution function is factorizable in longitudinal and 
transverse directions and that the shadowing can be neglected the average number of 
hard interaction 2xh(b, s) at the impact parameter b is given by : 

Xh(b, s ) = \ a jet{s) T N (b, s) (9) 

where T^(b, s) is the effective partonic overlap function of the nucleons at impact param- 
eter b. 

T N (b,s) = J d 2 b'p(b')p(\b-b'\) (10) 

with normalization / d 2 bTj^(b,s) = 1 and <jj et is the pQCD cross section of parton 
interaction or jet production J66[, |67j. Note that £ = 6/6 (s), where bo(s) provides 



a measure of the geometrical size of the nucleon 7t6q(s) = a s (s)/2 assuming the same 
geometrical distribution for both soft and hard overlap functions 

Xs(^s) = ^xo(0 (11) 

tt(e,*) = o^xo(e) (12) 

X^,s) = ^-[a s (s) + a jet (s)} Xo (0 (13) 
2a 

We note that x(£, s) is a function not only of £ but also of y/s because of the y/s 
dependence on the jet cross section <Jj e t(s) .Geometrical scaling implies on the other hand 
that Xs(£; s ) — Xo(0 is on ly a function of £ . Therefore, geometrical scaling is broken at 
high energies by the introduction of <x, ei (s) of jet production. 

The cross sections of nucleon - nucleon collisions can in this case be expressed as: 



CO 



a el = a (s) / d?(l- exp(- X (Z, s)) (14) 

Jo 

Pin = °v{s) de(l-exp(-2 X {i,s))) (15) 
Jo 



a tot = 2a (s) / deH-exp{- X ^s))) (16) 
Jo 

The calculation of these cross sections requires specifying a s (s) with a corresponding value 
of cut - off momenta p « 2 GeV/c flBSfl . 

In the energy range 10 GeV^ < \/s < 70 GeV, where only soft parton interactions are 
important, the soft cross section a s (s) is fixed by the data on total cross sections <y to t{s) 
directly. In and above the SppS energy range y/s > 200 GeV, a fixed cr s (s) = 57 mb and 
a mini-jet cutoff scale po = 2 GeV/c, leads to observed energy dependence of the cross 
sections and inclusive distributions. Between the two regions 70 GeV < y/s < 200 GeV, 
a smooth extrapolation for cr s (s) is used. 

In HIJING, a nucleus-nucleus collisions is decomposed into a sequence of binary col- 
lisions involving in general excited or wounded nucleons. Wounded nucleon are assumed 
to be q — qq string like configurations that decay on a slow time scale compared to the 
collision time of the nuclei. In the FRITIOF scheme wounded nucleon interactions follow 
the same excitation law as the original hadrons. In the DPM scheme subsequent collisions 
essentially differ from the first since they are assumed to involve sea partons instead of 
valence ones. The HIJING model adopts a hybrid scheme, iterating string-string collisions 
as in FRITIOF but utilizing DPM like distributions. In the SPS range the HIJING results 
for nuclear collisions are very similar to those of FRITIOF. However, HIJING provides 
an interpolation model between the nonperturbative beam jet fragmentation physics at 
intermediate CERN-SPS energies and perturbative QCD mini-jet physics at the highest 
collider energies ( RHIC , LHC ) . 

III. NUMERICAL RESULTS 

A. STRANGENESS IN PROTON - PROTON INTERACTION 

We used the program HIJING with default parameters: IHPR2(11)=1 gives the 
baryon production model with diquark-antidiquark pair production allowed, initial di- 
quark treated as unit; IHPR2(12)=1, decay of particle such as 7r° , , A , S , H , Q are 
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allowed ; IHPR2(17) = 1 - Gaussian distribution of transverse momentum of the sea quarks 
;IHPR2(8)=0 - jet production turned off for theoretical predictions denoted by HIJING 
model, and IHPR2(8)=10-the maximum number of jet production per nucleon-nucleon 
interaction for for theoretical predictions denoted by HIJING^ for comparison. 

In Table I the calculated average multiplicities of particle at Ei^ = 200 GeV in 
proton-proton (pp) interaction are compared to data. The theoretical values HIJING 
and HIJING {j) are obtained for 10 5 generated events and in a full phase space. The 
values HIJING {j) include the very small possibility of mini jet production at these low 
SPS energies. The experimental data are taken from Gazdzicki and Hansen |15] . 



The small kaon to pion ratio is due to the suppressed strangeness production basic 
to string fragmentation. Positive pions and kaons are more abundant than the negative 
ones due to charge conservation. We note that the integrated multiplicities for neutral 
strange particle < A >, < A >, < > are reproduced at the level of three standard 
deviations for pp interactions at 200 GeV . However the values for < p > and < A > 
are significantly over predicted by the model. This is important since as we shall see the 
A in S + S is significantly underestimated by HIJING. 

For completeness we include a comparison of hadron yields at collider energies ^/s = 
546 GeV( SppS -energies), for pp interactions, where mini-jet production plays a much 
more important role. From different collider experiments Alner et al.(UA5 Collaboration)) 
|86| attempted to piece together a picture of the composition of a typical soft event at the 
SppS |87j . The measurements were made in various different kinematic regions and have 
been extrapolated in the full transverse momentary ) and rapidity range for comparison 



as described in reference |86|. The experimental data are compared to theoretical values 
obtained with HIJING^ in Table II. It was stressed by Ward ]57[ that the data show 
a substantial excess of photons compared to the mean n + + 7r~ . It was suggested 
as a possible explanation of such enhancement a gluon Cerenkov radiation emission in 
hadronic collision pgfl . Our calculations rules out such hypothesis. Taking into account 
decay from resonances and direct gamma production, good agreement is found within 
the experimental errors. The experimental ratio ^ = 0.095 ± 0.009 is also reproduced 
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by HIJING^ model (0.099). We note that a detailed study of the ratios of invariant 
cross sections of kaons to that of pions as a function of transverse momenta in the central 
region was presented in [|67 . 



In the following plots the kinematic variable used to describe single particle properties 
are the transverse momentum p T and the rapidity y defined as usual as: 

1, E + p 3 E + p 3 

y=-ln- = ln (17) 

I h — ps rriT 

with E,p3 ,and tjit being energy, longitudinal momentum and transverse mass tjit = 



m o + Pt with m being the particle rest mass. 
In Fig. la, 3a, 4a, and 6a, we show rapidity and transverse momentum distributions for 
A 's (Fig.la,4a) and K^s (Fig.3a,6a) produced in pp scattering at 200 GeV. The theo- 
retical histograms obtained with HIJING (solid) and VENUS-4.13 (dashed) are compared 
with experimental data taken from Jaeger et al. The HIJING spectra for A , K® are 
close to the data at mid rapidity ||89|| , although the dip in the K® yield at mid-rapidity 
and the A peak in the fragmentation regions are not well reproduced (see also ref. flOPf). 
Unfortunately, more precise data are not available in pp interactions and those features 
could reflect experimental acceptance cuts. Similarly no detailed A spectra are as yet 
available in pp. 

In comparison with VENUS (taking 10 4 events) we note that this version seems to 
over-predict the pp — > A rapidity density at mid-rapidity by 50 — 100% in Fig. la, even 
though the rapidity integrated transverse momentum distribution in Fig 4a seems closer 
to the data. The yields in Figs. 3a and 6a are similar to those of HIJING with the 
dip structure in the data absent. 

The very sparse data base on pp strangeness production at SPS energies should be 
expanded in the future to improve the test of dynamical models before they are applied 
to the more complex nuclear collision case. Without A spectra in pp, for example, the 
need for the new dynamical mechanisms in that channel cannot be confirmed. 
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B. Multiplicities in pA and AA collisions 



In this section, we compare strange particle production in the HIJING and VENUS 
models to pA and AA data. Again we limit the study to A, A, K® to compare with recent 



data from Alber et al. [IS]. First we consider the average integrated multiplicities for 
negative hadrons < h~ > , negative pions < n~ > and neutral strange particles < K® > , 
< A > , < A > in pp,pS,pAg,pAu (minimum bias collisions) and SS, SAg, SAu (central 
collisions) at 200 AGeV. The default parameters of HIJING were used without mini-jet 
production (IHPR2(8)=0). The number of Monte Carlo generated events was 10 5 for 
HIJING and 10 4 for VENUS for pp,pA interactions and 5 • 10 3 for SS, and 10 3 for 
S + Ag, W, Au and PbPb collisions. 

The mean multiplicities are compared in Table III (for pp and pA interactions) and in 
Table IV (for AA interactions) with experimental data from Alber et al. ||18|| . Note that 
while the HIJING model describes well the integrated neutral strange particle multiplici- 
ties (except for < A >) in pp and pA interactions, there is a large discrepancy already for 
the light ion S + S reaction. 

It is worthwhile to mention that theoretical calculations have been done for pA 
'minimum bias' collisions and the experimental data are for the events with charged 
particle multiplicity greater than five, which contain a significant fraction (about 90 %) of 



the 'minimum bias' events [18 



In Table III and IV the data are compared also with other theoretical values obtained 
in some models : VENUS (as computed here), RQMD H, QGSM m, H and DPM 
models. The theoretical values DPM 1 are from the Mohring et al. f49fl , version of DPM 
which include additionally (qq) — (qq) production from the sea into the chain formation 
process and the values DPM 2 are from the Mohring et al. ||50|| , version of DPM which 
include chain fusion , as a mechanism to explain the anomalous antihyperon production. 

Alber et al. [18[] have considered that the total production of strangeness should be 
treated in a model independent way using all available experimental information for ratio 
Es expressed as : 
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< A > +4 < K° s > 

b 3 < 7T > 



(18) 



We have calculated this ratio in HIJING approach for the above interactions and the 
corresponding numerical predictions are shown in Table V. We note that there is much 
less discrepancy between HIJING and the data for this particular ratio. We conclude 
from this that such ratio is insensitive to the underlying physics and therefore should 
NOT be used for any further tests of models! This ratio hides very effectively the gross 
deficiencies of the HIJING model in SS reactions pointed out later in the comparison to 
the rapidity and transverse momentum distributions. We include Table V only to prove 
the futility of studying the systematics of such ratios in the search for novel dynamics in 
nuclear collisions! 



C. Single inclusive distributions for neutral strange particles in pA and AA 

The main results of the present study, are contained in Figs. 1-6. Figure 1 is our most 
important result revealing the systematics of A enhancement from (a) pp to (d) SAu. In 
part (a) the pp data at mid rapidity are seen to be well reproduced by HIJING. However, 
the new minimum bias pS data |H| in Fig. lb clearly shows a factor of 2 — 3 discrepancy 
with respect to the linear extrapolation from pp as performed by HIJING. The effect of 
double string fragmentation and final state cascading, as modeled with VENUS is seen 
on the other hand to account for the observed A enhancement. We note, however, that 
in pp, VENUS over-predicts the A yield at mid rapidity. Some fraction of the aggreement 
in pS with VENUS may be due to this effect. The overprediction of midrapidity A's in 
pp by VENUS was shown in Fig. 10.20b of ref. ||60|| , but was not emphasized there. If 



both the pp and pS data on A production are correct, then the most striking increase of 
hyperon production therefore occurs between pp to pS reactions. 

The strangeness enhancement in minimum bias p + S is striking because the number 
of target nucleons struck by the incident proton is on the average only two! The step from 
single p+p to triple p+p+p reactions therefore apparently leads a substantial enhancement 
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of midrapidity A's which obviously cannot have anything to do with equilibrium physics. 

In central S + S reactions shown in Fig. lc, the discrepancy relative to HIJING grows 
by another factor of two. We note that the new data ||18|| shown here have increased 
substatially relative to earlier data |7|], || due to inclusion of lower transverse momentum 
regions and A's originating from the decay of E, S in the analysis. Including these decay 
channels, VENUS is seen to reproduce the new data as well. We note that with RQMD 
the excess A's is also reproduced with the introduction of rope formation (see Table IV). 
For heavier targets, Ag, Au , in Fig. Id, the discrepancy relative to HIJING is in fact less 
dramatic than in S + S . 

In central S + S , on the average each projectile nucleon interacts with only two target 
one, but each target nucleon also interacts with two projectile ones. In effect, then S + S 
reactions probe strangeness production in four nucleon interactions p + p + p + p — > 
A + X. Such reactions appear to be approximately four times as efficient in producing 
midrapidity A's as two nucleon interactions in part (a). Our main conclusion therefore 
is that strangeness enhancement is a nonequilibrium dynamical effect clearly revealed in 
the lightest ion interactions. 

Further support for this conclusion is shown in Figs. 4 a,b,c, where the transverse 
momentum distributions are compared. We see that there is an enhancement of the A 
transverse momentum relative to pp in pS. Comparing to VENUS we can interpret Fig. 
4b as evidence that the enhanced transverse momentum of A in pS is due cascading. The 
discrepency in Fig. 4c between VENUS and the data in SS may be due to the rapidity 
cuts in the data, which we have not included in the calculated spectra. In all cases the 
deficiency of the linear extrapolation via the HIJING model is clearly evident. For heavier 
targets, S + Ag, Au, the transverse momentum distribution predicted by VENUS is close 
to the data. 

The same general conclusion emerges from the systematics of A and K® production s 
in Figs 2,3 and in Figs 5,6 respectively. In Fig. 2a, the agreement between HIJING and 
VENUS and the data on the p + S — > A must be viewed with caution since as shown in 
Table III, both models overpredict the integrated A multiplicity by a factor 2 — 3. Given 
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the absence of more detailed rapidity and transverse momentum distributions for A, it 
is not possible to determine whether the pp and pA data are compatible. However, at 
least the step form pS to SS in Fig 2b indicates a possible factor of two enhancement of 
A similar to the comparison of Fig. lb,lc for A. As in the case of A production, there 
appears to be no further A enhancement from SS to SAu. As regards to the transverse 
momentum distributions in Fig. 5, we note that as in Fig. 4 the A emerge with higher 
p± in pS than in pp in accord with the VENUS model. We note that in Figs 5 c,d, the 
norm theoretical curves is obtained intergrating over the full rapidity interval, while the 
norm data are limited to a smaller domain as shown in Fig. 2 b,c. 

In the case of K® production in Figs 3, 6, the same general trends are seen but in a 
less dramatic form. 

We conclude that the new data indicate that the origin of strangeness enhancement 
in heavy ion collisions may be traced back to non- conventional and necessarily non- 
equilibrium dynamical effects that arise in collisions of three or more nucleons. However, 
this conclusion is forced upon us by the sytematics of the new light ion data on p + S and 
S+S reported in |18| . As shown in Figs lb, 2a, and 3b, those systematics, especially in pA, 
differ considerably from the trends of earlier NA5 data [TJJ] and preliminary NA36 data 



liTj - f26| . Those data for heavier target nuclei incidate substantially less enhancement of 
midrapidity A, A, K® than do the NA35 data on p + S. Part of the difference between 
these data sets may be due to different acceptance cuts and the inclusion or rejection of 
fragments from decay of higher mass hyperons. Obviously, the difference between these 
data sets must be resolved. Until then, the NA36 data must be regarded as an important 
caveat on our conclusions. 

For completeness we show also in Figs 7 the linear extrapolations of HIJING to Pb+Pb 
at 170 AGeV for all positives (Fig 7a) and all negatives charges (Fig 7b), for A (Fig 7c) 
and for A (Fig 7d). It will be interesting to compare these extrapolations with upcoming 
data to test if the strangeness enhancement increases from SS to PbPb. 

We include the two dimensional distributions in Fig. 8 to emphasize that strangeness 
enhancement analyses restricted to narrow rapidity and transverse momentum cuts, es- 
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pecially with simplistic fireball models, may completely miss the global non-equilibrium 
character of the data. 



IV. CONCLUSIONS 

In this paper we performed a systematic analysis of strange particle production in 
pp ,pA and AA collisions at SPS CERN-energies using the HIJING and VENUS models. 
The most surprising result is that the breakdown of the linear extrapolation from pp 
data to nucleus-nucleus in the strangeness channel already occurs in minimum bias pS ! 
The apparent enhancement of A, A and K® at midrapidities in pS reactions by a factor 
of 2 indicates that the mechanism for strangeness enhancement in heavier ion collisions 
must be associated with non-equilibrium dynamics involving multiparticle production and 
not with equilibrium quark-gluon fireball. In minimum bias pS one projectile nucleon 
interacts on the average with only two target ones. The data on pS therefore indicate 
the existence of new dynamical mechanisms for strangeness production that becomes 
operative in p + p+p collisions. The new data |18j on central S + S show another factor of 
2 enhancement of strangeness production relative to pS. This light ion reaction basically 
probes multiparticle production in p + p + p + p. The strangeness enhancement in heavier 
target systems apparently saturates at the S+S level. We also showed that traditional 
analysis of strangeness enhancement in terms of ratios of integrated multiplicities is very 
ineffective since those ratios hide well defects of the detailed rapidity and transverse 
momentum distributions predicted by models. 

The agreement with VENUS and RQMD results suggests color rope formation as a 
possible mechanism. However, to clarify the new physics much better quality data on 
elementary p + p as well as on other light ion p + a,C, S and a + a, C, S reactions will be 
needed. Especially, the discrepancy between NA35 and NA36 must be resolved. Only then 
can strangeness enhancement systematics used meaningfully in the search for signatures 
of quark-gluon plasma formation in future experiments with Au + Au and Pb + Pb. 
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FIG. 1. Rapidity distributions of A produced in pp interactions at 200 GeV (Fig. la). The 
data for pp (black small circles) are from Jaeger et al. [p9|| . Rapidity distributions of A produced 
in minimum bias pS (Fig. lb) and central 55 (Fig. lc), 5 Ag (Fig. Id) and SAu (Fig. Id) collisions 
at 200 AGeV. HIJING and VENUS results are shown by solid and dashed histograms (for 
pp ,pS , 55 , SAu). The new NA35 data (pS , SS -full circles; SAg - stars ; SAu -full triangle) 



are from Alber et al. [18]. The open circles show the distributions for SS collisions reflected at 



Uiab = 3.0 . In Fig lb, earlier NA5 data onp + Ar (open diamond) from ref. [14| and preliminary 
data on p + Pb (open squares) from NA36 fll9|p6f are shown for comparison. In Figlc NA36 
data on 5 + 5(open squares) are also shown for comparison. 

FIG. 2. Rapidity distributions of A produced in minimum bias pS (Fig. 2a) and central 
55 (Fig.2b), 5^ (Fig.2c) and SAu (Fig.2d) collisions at 200 AGeV. Solid and dashed his- 
tograms are as in Fig.l . The NA35 data (full circles) are from Alber et al. [18|. The open 



circles show the distributions for SS collisions reflected at yi a b = 3.0 . In Fig. 2a, the open 
squares correspond to preliminary p + Pb data from NA36 [19|,26]. 



FIG. 3. As in Fig.l but for K° s particles. 

FIG. 4. Transverse momentum distributions of A produced in pp interactions at 200 AGeV 
(Fig. 4a). The data (black small circles)for pp interactions are from Jaeger et al. |89|| . Trans- 
verse kinetic energy distributions of A produced in minimum bias pS (Fig. 4b) and central 
55 (Fig.4c), 5^5 (Fig.4d) and SAu (Fig.4d) collisions at 200 AGeV. HIJING and VENUS re- 
sults are shown by solid and dashed histograms resp. (for pp ,pS , 55 , SAu). The NA35 data 
(pS , 55 -full circles; SAg - stars ; SAu -full triangle) are from Alber et al. [18|. 



FIG. 5. Transverse kinetic energy distributions for A particles produced in minimum bias 
pS (Fig.5a) and central 55 (Fig.5b), SAg (Fig. 5c) and SAu (Fig.5d) collisions at 200 GeV per 
nucleon. Solid and dashed histograms are as in Fig. 4. The experimental data (full circles) are 



from Alber et al. |18|. 
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FIG. 6. As in Fig.4, but for K° s particles . 

FIG. 7. Predicted rapidity distributions for central (5 = — 1 fm) PbPb collisions at 170 
AGeV with the HIJING model for all positive charges (Fig.7a), all negative charges (Fig.7b), 
A (Fig. 7c) and A (Fig.7d). 

FIG. 8. Unnormalized rapidity y and transverse momentum px distributions for A 
(Fig.8a,b) and A (Fig.8c,d) for central PbPb at 170 AGeV from HIJING. 
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TABLE I. Particle multiplicities for pp interaction at 200 GeV are compared with data from 
Gazdzicki and Hansen 



PP 
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TABLE II. Particle composition of p + p interactions at 540 GeV in cm. 



Particle type 


< n > 


Exp. data 
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TABLE III. Average multiplicities for negative charged hadrons and neutral strange hadrons 



in pp and pA interactions. HIJING and VENUS model results are compared with others recent 
estimates using RQMD, QGSM, DPM and with data from Alber et al. @. 



Reaction 
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TABLE IV. Average multiplicities for negative charged hadrons and neutral strange hadrons 
in A A interactions. HIJING and VENUS model results are compared with others recent esti- 



mates using RQMD, QGSM, DPM and with data from Alber et al. 



Reaction 
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TABLE V. The mean multiplicities of negative pions and Es ratios(see the text for defini- 
tion) for nuclear collisions at 200 AGeV. The data are from Alber et al. [jl8|and the NN data 



are from Gazdzicki and Hansen 15]. 
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P + P,P+S(200 GeV),S + S,S+Au(200 AGeV) HIJING(-) vs. VENUS( — ) 
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